Generating quantizing pseudomagnetic fields by bending graphene ribbons 
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We analyze the mechanical deformations that are required to create uniform pseudomagnetic 
fields in graphene. It is shown that, if a ribbon is bent in-plane into a circular arc, this can lead to 
fields exceeding lOT, which is sufficient for the observation of pseudo-Landau quantization. The arc 
geometry is simpler than those suggested previously and, in our opinion, has much better chances 
to be realized experimentally soon. The effects of a scalar potential induced by dilatation in this 
geometry is shown to be negligible. 



Graphene exhibits a number of unique features not 
found in conventional metals and insulators.^ ^' Among 
them is the possibility to stretch graphene elastically 
by more than 15 %'^, and to co ntrol in different ways 
the induced strainPSsEISllIinilll. This offers a prospect 
of tuning electronic characteristics of graphene devices 
not only by external electric field but also by me- 
chanical strairPEH^ a p ossibility being extensively dis- 
cussed theoreticalljff^3 i4 15 16 17 is 19 20 particu- 
lar, the presence of two valleys at the opposite corners of 
graphene's Brillouin zone implies that long wavelength 
lattice deformations induce an effective gauge field acting 
on the electrons and ho les, which has the opposite sign for 
the two valleys . ^^^ l ^^ l This yields an enticing possibility 
of creating such gauge fields that would mimic a uniform 
magnetic field B and, consequently, generate energy gaps 
in the electronic spectrum and lead to a zero- 5 analogue 
of the quantum ffall effect. Both isotropic and uniaxial 
strains resultP^'^il jjj ^ero pseudomagnetic field Bs but 
as shown recentlj^ deformations with a triangular sym- 
metry can lead to strong uniform Bg. Moreover, the 
strained-induced pseudomagnetic field can easily reach 
quantizing values, exceeding lOT in submicron devices 
for deformations less than 10%.'^ 

Unfortunately, all the geometries of applied strain, 
which were suggested previously,'^ are rather difficult to 
realize experimentally. In this Communication, we report 
an alternative strain configuration that does not require 
a complex triangular symmetry and, in fact, is a straight- 
forward extension of the geometry ty pically used in ex- 
perimental studies of strained deviceaSU^UJJ. We have 
found that simple in-plane bending of graphene ribbons 
(see Fig.[T]) should lead to strong practically uniform ^5. 
We believe that this finding can speed up the observation 
of the pseudomagnetic quantum Hall effect and related 
phenomena. 

First, let us complete the analysis of Ref.24 by clas- 
sifying the strain distributions that are compatible with 
equilibrium elasticity and lead to a uniform pseudomag- 
netic field. 

We will use the coordinates that are fixed with respect 
to graphene's honeycomb lattice in such a way that the 




FIG. 1: (Color online). Sketch of the suggested bending ge- 
ometry that would generate a uniform pseudomagnetic field 
and open band gaps in graphene's electronic spectrum. The 
graphene rectangle (lower image) is bent into a circular arc 
(upper). 



X axis corresponds to a zigzag direction. In this case, the 
gauge field A acti ng on charge carriers in graphene can 
be written a^^IIll] 

Ax — ic {uxx '^yy) 

Ay = ^2c^Uxy (f) 

where (3 = ~dlog{t) /dlog{a), where t « 3eV is the 
electron hopping between orbitals located at nearest 
neighbor atoms, a « I.4A is the distance between them, 
c is a numerical constant that depends on the details of 
atomic displacements within the lattice unit cell, and Uij 
is the strain tensor. The two signs correspond to the two 
valleys, K and K' in the Brillouin zone of graphene. 

In two dimensional elasticity problems, it is convenient 
to study the stress tensor, atj — dT /duij, where T is the 
elastic energjI^H. The gauge field can be written in terms 
of the stress tensor as 

Ax — ic^^^ {^^xx ^yy^ 

Ay — Cxy (2) 
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where ^ is a Lame coefficient. Furthermore, possible 
stress distributions that describe two dimensional elastic 
systems in equilibrium can be written in terms of com- 
plex variables z = x + iy and z — x — iy ai 

_ d^f{z,z) 
dx-^ 

d'f{z,z) 
dxdy 



'vv 



(3) 



Here /(z, z) is either the real or the imaginary part of a 
function 



J^{z, z) = J^i{z) + zJ^2{z) 



where 3'i{z) and T2{z) are analytic functions. For the 
case of pure shear deformations considered in ReflMI 
^2 = but here we will not restrict ourselves by this 
limitation. Since both stress and A are given by the sec- 
ond derivatives of T , whereas Bs is given by the first 
derivatives of ^, a uniform Bs necessitates T to have a 
cubic dependence on the coordinates. Such a function 
must have the following structure: 

T{x, y) ^ci{x + iyY + C2{x - iy){x + lyY (5) 

where ci and C2 are arbitrary constants. Separating the 
real and imaginary part of equation (6), we find four 
possible functions that result in uniform Bg 



x'^ — 3xy^ 
x^y + y^ 



(6) 



The second pair of the solutions is equivalent to the first 
one by swapping the axes. For the lattice orientation 
used in Eq. ([T]), the first pair leads to a gauge field such 
that Ax (X X and Ay cx y and, accordingly, Bs is zero. 
Hence, the stress distributions that give rise to a uni- 
form pseudomagnetic field can be expressed in terms of 
a superposition of the functions in lines 3 and 4. The 
strain configuration found in rcf 24 involves only the 3rd 
function f{x,y) cx 3x^y — y^, which leads to a unique 
solution for the shape of graphene fiake where such dis- 
tribution of stresses can be created by normal forces only. 
Unfortunately, this solution is not easy to realize exper- 
imentally. The use of both 3rd and 4th functions offers 
further possibilities. 

In the following, we consider the deformations required 
to create a uniform Bs inside a rectangular graphene 
crystal, of width W and length L, with normal forces 
applied at the left and right boundaries as sketched 
m Figure [2| The non-deformed crystal fills the region 
-L/2 < X < L/2,-W/2 < y < W/2. Let us write the 
forces at the right and left boundaries x = zLL/2 as 



R,L 



fa + fi y 



R.D 



(7) 



FIG. 2: (Color online). Stretching geometry leading to a 
uniform pseudomagnetic field inside a rectangular graphene 
sample. Normal forces are applied at two opposite boundaries 
and their magnitude is indicated by the length of the plotted 



The condition of zero total force and zero total torque 
requires = -/(f = /q and //^ = = /i where /o 
and /i are constants. The absence of forces at the upper 

there. At 



and lower edges implies that ayy 
the right and left edges, we have 



o-xx ^ [ y + 

CJxy = 



/o 



(8) 



where ctq is a constant that depends on applied forces. A 
stress distribution within the crystal, which is compat- 
ible with these boundary conditions, is generated by a 
function f{y) = Ciy^S + (/o//i)(2/V2)]. This function 
can be considered as a superposition of solutions 3 and 
4 of Eq. ([6]), which leads to a uniform Bs, and a con- 
stant term that describes a uniaxial strain and does not 
give rise to an additional pseudomagnetic field. The lat- 
ter term ensures that the lattice is stretched everywhere, 
there is no possibility for out of plane deformations. In- 
side the rectangular crystal, axx depends only on y in the 
manner described by eq. ([8|, and ayy = Uxy = 0. From 
this stress distribution, we find the lattice distortions 



Urr. = uq ( 2xy + j-x 



A 



Wo 



A + 4/1 



y 



fo 



■y 



(9) 



where uq is a constant that defines the maximum stress. 
These displacements lead to the curved shape shown in 
Fig. [2j which was drawn using the reportecr^ Lame co- 
efficients of graphene, A ~ 3.3eV and /i w 9.4eV 
A-2. The maximum strain occurs at the top and bottom 
boundaries and can be estimated as Umax ~ uo(W + 
fo/fi)- The pseudomagnetic field inside the graphene 
crystal is given by 



2^>oUo 

Bs = c/3 = C/3-— 

a aW 



(10) 



where $0 is the fiux quantum. The effective magnetic 
length is £b = \/ (aW) / (/3w). This field has the same 
dependence on the crystal dimensions and the maximum 
strain as in the examples discussed in Ref,24,. For W « 
0.1 micron and u « 10% the generated effective field is 
of the order of 20T. 
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Experimentally, it may be difficult to create the precise 
stress distribution prescribed by Eq. ([s]). However, one 
can see that the required shape of the strained graphene 
crystal in Fig. 2 resembles an arc of a circle. To this end, 
we consider next the geometry in which a rectangular 
graphene crystal is bent into a circular arc, as sketched 
in Fig. 1 and shown in more detail in Fig. [3^. Note that 
this geometry is in fact standard for experimental studies 
of strain (see, e.g., Ref ]8ll0irT|l with the only difference 
that the bending should be applied in plane rather than 
out of plane of a graphene sheet. 

If the radius of the inner circle defining the lower edge 
in Fig. 3a is R, the shape of the deformed rectangle is 
given by 



Ux{x,y) = (i? + y)sin 
Uy{x,y) = (R + y) cos 
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The undistorted rectangular shape is recovered for 
L/R 0. The displacements in Eq. (11) can be ex- 
panded in powers of R~ 



and the leading terms are 
Ux{x,y) 



xy 
R 



Uy{x,y) 



X 

2R 



(12) 



These displacements do not exceed max(i, W)'^/R. The 
next terms lead to corrections bound by max(i, W)^ /R^. 
The distortions in Eq. ( [T2| coincide with those in eq. ([9| 
in the limit of vanishing Poisson ratio A/(A + 2^) — > 
and lead to a uniform Bs inside the sample. The maxi- 
mum strain is L/2R. For an arbitrary value of L/R, the 
pseudomagnetic field is given by 



Bs{x,y) = -4c 



aL 
R + 



■ ( 




'2x 






arcsm — 


cos 


— arcsm 




X 


\2Rj 











arcsm 



L 
2R 



(13) 



to Bs ~ 
( 10 ) and the es- 



For L/R — > 0, the field reduces 
— c(/3<[>o)/(ai?), in agreement with Eq. 
timates given in Refl24l The relative corrections to the 
constant value of Bs are of the order L/2R, that is, the 
maximum strain. An example of the field distribution 
described by Eq. ( 13 1 is plotted in Fig. [sb. 



The found strain is not purely shear but it also contains 
a dilatation. The latter gives rise to an effective scalar 
potentiaP^, in addition to the discussed pseudomagnetic 
field. Below, we show that due to screening, the extra 
potential does not radically affect Landau quantization. 

Following RefHU and using eqs. [9] this potential is: 



V{x, y) = Vo {d,j,Ux + dyUy) = Vq 



2u„y 



/i(A + 4m) 

(14) 
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FIG. 3: (Color online), (a) - Rectangular graphene sample 
deformed into an arc. The radii of the lower and upper edges 
are R and R+W, respectively. The plot is for R — 5 x L. (b) 
- Effective magnetic field, in Teslas, for the same deformed 
ribbon. Dimensions are W = 200nm, L — 192nm, R = 5 x 
L — 960 nm. The maximum strain is 10%. 
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FIG. 4: (Color online). Variation in density, 5p(x,y), in di- 
mensionless units (see text) due to the screened scalar poten- 
tial induced by strains. The aspect ratio is L/W = 1. 



where Vq « 3eV, estimated from the linear rise in the 
work function of graphene under compression between 
and 10% straiifi^ (note that in Refl^a much larger value 
of 16 eV is quoted, from old experimental data on trans- 
port properties of graphite) . The constant contribution 
gives a rigid shift to all the levels while the non uniform 
term is equivalent to an effective electric field along the 
y direction. Unlike the strain induced gauge field, this 
potential will induce a charge redistribution and will be 
screened by the carriers in graphene. 

We consider first the screening expected if we assume 
that the flake is a perfect metal, neglecting corrections 
due to quantum properties of the electron gas. The in- 
duced charge density, 5p{r) can be thus obtained from 
the condition 



Spir') 2Vqu 



W 



5p{r)y 



(15) 



This equation can be rescaled by the substitution 
5p{x/L,y/W) = {W+L)/{2L)x{2Vau)/{e^W)x5p{x,y). 
where the function 5p{x, y) depends only on the aspect 
ratio, W/L. 

The density of carriers in a given Landau level due 
to the effective field given in Eq. [lO] setting c = 1, is 
PLL = 2{aW)/{Pu) so that 6p/pLL ^ {W + L)/{2L) x 
(4Voa)/(/3e^) x 6p{x/L,y/W). Hence, the non uniform 
density induced by the effective field, in units of the 
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difference in densities between different quantum Hall 
plateaus, is approximately given by 5p shown in Fig. [4] 

To provide an ideal metallic screening, the chemical 
potential should coincide at each point with the posi- 
tions of one of the Landau levels. Hence, the quantum 
energy of the carriers is larger than in a uniform electron 
distribution, when the Fermi energy lies in a pseudogap 
between the Landau levels. Thus, in addition to the clas- 
sical screening energy, discussed above, we also must add 
the change in energy due to the changes in occupancies 
of the electron levels in the presence of the scalar poten- 
tial. We first assume that the induced electronic density 



is given by the solution of Eq. 15 and that the electronic 
states are the Landau levels induced by the effective field 
in Eq. |10| Then, the two contributions to the energy, for 
W ^ L, are of the order: 



quantum 



mi) 



(16) 



The calculations shown in Fig. |4] suggest that {\Sp\) = 
f « 10~^ - 10"\ and ~ ~ Voa. Then, the scalar 
potential is screened, and the process can be de scribed 
by the classical model outlined earlier, if y/L/a > /"^, 
that is, L > 10^ — lO'^nm. For smaller sizes, the rigidity 



of the quantum levels induced by the gauge potential 
prevents any rearrangement of the charge inside the flake. 
A detailed theory of screening in this situation will be 
presented elsewhere. In either case, the electronic states 
are well described by the effective Landau levels induced 
by the field in Eq. 

To create the required strain experimentally, one can 
think of depositing graphene ribbons onto a rectangular 
elastic substrate and deform it in the manner prescribed 
by eq. (9) or by bending it into a circular arc (Fig. 1). 
Crystals rigidly attached to the substrate can then be 
of arbitrary shape, as the strain distribution in the sub- 
strate would project onto graphene and give rise to a 
(nearly) uniform Bs- Note that macroscopic substrates 
would require the use of rubber-like materials capable of 
withstanding very large strains, such that local deforma- 
tions projected on a submicron graphene crystals could 
still reach w 10%. 
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